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tfie present invention reiat09 to a method for locatins a eased bore* 
hole from an adjacent open or uncased borehole • Many tines during the drill- 
ing of oil wells it is necessary or desirable to locate a previously drilled 
well, i&ich is cased with a ferromagnetic material » from a second well which 
is ancased. For example, in the case of a blowout in a previously drilled 
deep well, it nay be more economical to intercept the previous well at a 
shallow depth above the total depth and kill the formation than to drill to 
the total d^th of the original well. Likewise* when drilling multiple wells 
from a single offshore platform, it is desirable to know the exact location 
10 of the existing wells with respect to the well that is being drilled to |»re- 
vent an accidental interception of previously drilled wells. 

In the past it has been impossible to accurately predict the loca- 
tion of a previously drilled well ftom a second well. Hhile the need for such 
a device has long been felt in the petroleum industry, it has not been satis- 
fied by available equipment. 

Object of the invention is a mdthod for detexnining the distance 
and direction to a cased borehole from measurements made in adjacent uncased 
borehole, said cased borehole being cased with a material having a remanent 
magnetization. 

20 According to the invention there is provided a method for deter- 

mining the distance and direction to a cased borehole from measurements made 
in adjacent ancased bordiole, said cased borehole being cased with a material 
having a remanent magnetization, said method comprising: measuring, in the 
uncased borehole at a plurality of vertically spaced positions « the total 
strength of the magnetic field existing at each of the positions; measuring, 
in the uncased borehole at substantially the same positions, the deflection 
of a magnetic compass from o geographical reference; measuring the depth of 
the spaced positions from a known reference point; assuming a value and 
direction for the earths magnetic field at the positions; and determining 

SO by machine calculations the distance and direction to the cased borehole 

using said measurements and assumed value for the earth's field and location 
for the cased borehole, said calculations using a least squares fit of the 
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measured values to calculated values in the fomula for the nagnetic Induc- 
tion surrounding a long cylinder. 

Ihe present invention provides a method for locating a previous ly 
drilled well whidi is eased with a material having a remanent magnetisation, 
to have a remanent magnetisation in the casing is not difficult since most 
well casing is electrooagnetlcally inspected before it is Installed. The 
eiectranagnetic inspection leaves a remanent magnetisation in the casing. 

Since casing is normally installed in individual sections that are 
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joined together, the remaoent oagnetlzation of unperturbed eaeing nonially is 
periodic. 

Magnetonatera are avial^la for meaauring the total strength of a 
magnetic field. The field measured hy the magnetaneter vllXy of course, be 
the eoDbinfttion of the field trm the aagnetised casing plus the field of the 
earth, Xnovlng the total joagnetie field and the earth's field, one can cal- 
culate possible locations of the previously drilled iiell by t^tmAr.^ various 
locations and ccnputing the JBinimnm least squares error of fit to observed 
total field data. After wilmlfffting possible locations of the previously cas- 
10 ed well and knoviog the etrengtb of the earth's field, one can determine the 
strength of the field contributed by the cased nell. Knowing the location 
and strength of the field contribnted by the casing, one can calculate the 
dflfleetion of a nagnetic conpass caused by the cased vell« Xhe calculated 
deflection can then be coos^ared to tb» aetual deflection of the inagnetic coo- 
pass as datexsdned by a survey that includes both a gyro and magnetic coopass. 
Fran a eonparison of the oaleulated deflection Tilth the actual deflection, one 
can then datendne the location of the eased neU. In those cases nhere one 
is only interested in taowing the location of the first well so thet the second 
weU does not intercept it, one can utiHae the infomation to control the 
80 direction of the second well. 

The present invention will be ctescribea by way of exan^le in more 
detail with reference to the attadiea Orawinga in which: 

Fl«;ttre 1 is an isometric view of the locations of the two vella; 

Figure 2 is a representative plot of the total inagnetic fiela measur- 
ed by a magnetometer, the magnetic and gyro compass headings; ana 

Figures 3-5 are representative plots of the aefljoctiQn of a magnetic 
coD^ass and calculated defleotiona based on the measured magnetic field. 

The theory of magnetic fields including both the earth' a f ieln and 
an additional field is well understood by those skilled in physics. This theory 
30 can be ^iied to the case of a hollow, infinitely long, magnetised cylinder 
extending into the earth. When tha megnetlaea cylinder is a well casing , one 
can use the theory or magnetic fields to derive the following fonnulaa for the 
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. h (2nirp/L) cos (2nvZ/l) ♦ b h (2mp/l) sin (2niiZ/L) 



magnetic field due to a magnetised well casing: 

(1) B a 

n«i 

^ K K 

(2) B- =^ a^ ^)(2m^p/L3 sin {2niiZ/L) - b '^of2nirp/» cos (2nii2/L 

n»i " lC^(2niia/L) K^(2nifa/i) 

In the foroala^ is the radial conqponent of the external field 
and B^ is axial conponentj a is the outer radius of the casing and L is the 
wave length of the periodic magnetisation of the casing* and are modi- 
fied Bessel Atnctions of the second kind, a and b are constants^ p is the 

n n 

distance from the axis of the easing to the point of measurement, and Z is 
the vertical distance. 
10 Formulas 1 and 2 above apply to the case where the magnetised casing 

has a regular periodicity in the magnetization*. Normally* the casing vill be 
subjected to only a single stress since the casing is suspended from a hanger 
member in the borehole and remains in tension. In the case of a blowout or 
other abnormality, the casing can be subject to factors that distort the per- 
iodic character of the magnetic field external to the casing. In this- case, 
L is chosen larger than the interval over which data is being fit and fon- 
ulas 1 and 2 can still be used. 

In Figure 1, there is shown a schematic three-dimensional view of 
the two wells, one having a magnetized liner 10 and the other well 11 being 
20 drilled to intercept the first. Ihe liner 10 is ferromagnetic and has a re- 
manent magnetization* The remanent magnetization may result from an electro- 
magnetic inspection of the casing before it was installed or may be purposely 
provided by other means. The horizontal coordinate system assumes that the 
first well is located at the axis of the coordinate system and the second well 
is located to the east of the first well and drilled at an angle to inter- 
cept the first well. At known depths in the second well, it is assumed that 
the first well can be in the various positions along a line 12, such as IS, 
while the second well is in the position 18. Thus the distances from the 
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assuoed positions of tho first woll to the second well are H and 0 or they 
can be represented by the single line p which is the shortest distance from 
the axis of the second well to the assumed position of the first well in the 
horizontal plane. The north-south and east-west aagnetic axes are also 
shown. 

Tlie method of the present invention consists of the steps of first 
measuring in the second well the total magnetic field. As pointed out above, 
this measurenent can best be made using a sensitive magnetometer; for example, 
a magnetometer having a sensitivity of at least lO"^ gauss would be suitable. 
At the same time that the magnetic field is measured » the particular depths in 
the second well should also be measured* In addition, the deflection of a 
magnetic compass in relation to a gyro compass or similar device should also 
be measured. 

From the total measured field, one can proceed to calculate the 
strength of the radial field and vertical field B^ resulting from the mag- 
netized casing of the first well at the position of the second well. The 
earth's field is assumed to have a value of about .55 gauss and a direction of 
spproximtoly 65** from the horizontal and directed along an azimuth corres- 
ponding to magnetic north. The exact values for the earth's field will vary 
at various geographical locations. Those given above are approximate values 
f6r the central Mississippi area of the United States. The total measured 
field can be expressed vectorially as follows: 

Wherein B is the total field, B„ the earth's horizontal field, B^ the magnet- 
ized casing's radial field; B^ the earth's vertical field; B^ the magnetized 
casing's vertical field; a the angle between magnetic North and a line join- 
ing the two wells; and e^, e^ and e^ are unit vectors in the Iforth, Bast and 
vertically downward directions, respectively. 

If both sides of the above vector equation are squared, one obtains: 
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Prom an assuaed position fbr the firat woll, values for p and a can 

bo obtained and the unknowns a and b in the foraulae 1 and 2 calculated . 

n n 

This calculation is dono with L greater than the length of the path over which 
the calculations are being considered. Also n should be varied to provide 
sufficient wave lengths to determine the best fit of the observed data to the 
calculated data. The above calculation is carried out by first considering 
only the first order terms , 2 B^B^ Cos a ♦ 2 B^B^ to obtain preliminary value 
of the unlcnowns a^ and b^. These values ore then used to calculate the values 
of the s<iuared terms ♦ B^^] ♦ [b^^ ^ ^® calculated 

values of the squared terms, new calculations are made for the first order 
terms. The calctdations are continued until the best fit of the calculated 
value of the total field and the observed value of the total field is obtain- 
ed. Hie best fit can be obtained using a conventional least squares fit of 
the measured data to calculated values, i.e., B^ measured -B^ calculated. 

While no methods have been described in detail for carrying out the 
above calculations, anyone skilled in the eoi^ter art can program a computer 
to provide a least squares fit of measured data to calculated data. The use 
of least squares to fit measured data to calculated data Is described in many 
books on applied mathematics. Hius, no description of a program for a parti- 
cular computer is included. 

After the strengths of the radial and vertical fields in the second 
well resulting from the magnetised liner in the first well have been calcul- 
ated, the resulting deflection of a magnetic compass can be calculated. Ihe 
calculated values of deflection can then be compared to the measured values 
of deflection. This is important since positions of the first well on 
opposite sides of the second well may give substantially the same errors, 
ifhile the least mean square errors are the same, the deflection of magnetic 
compass will be materially different at the two positions. Shown in Figure 3 
is the apparent magnetic declination derived from measurements of the azimuth 
of the second well with magnetic and gyromultishot survey instruments. The 
data are representative of that obtained in a well . Figures 4 and 5 represent 
magnetic declinations calculated for two assumed locations of the first well. 
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It can be seen that the assuned locotion for the calculated declinations 
shown in Figure 4 is the location of the first well. 

The ahove procedure was used to intercept a well which had been 
drilled to approximately 20,000 feet and cased with a laaterial having 
renaneat nagnetization. Th& well subsequently blew out and it was desired 
to intersect the well at approxiaately 10,000 feet to avoid the drilling of 
a relief well to 20,000 feet to kill the formation. Also, it was believed 
that the casing in the original well had parted at a few thousand feet and 
placed the casing at the depth of interest in compression. The data shown 
in Figure 2 represent the magnetometer data obtained over an interval of 
approximately 500 feet in the second well using a magnetometer of high 
sensitivity. 
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Tim EMBODIMEhrrS OF THE INVTENTION IN WHICH AN EXCUISIVE 
PROPERTY OR PRIVILEGE IS ClAIMFiD ARB DEFINED AS FOLLOWS: 



!• A nethod for determining the distance and direction to a cased bore- 

hole fipon measurements made in adjacent uncased borehole, said cased borehole 
being cased with a material having a remanent magnetisation, said method com- 
prising: measuring, in the uncased borehole at a plurality of vertically 
spaced positions, the total strength of the magnetic field existing at each 
of the positions; measuring, in the uncased borehole at substantially the 
same positions, the deflection of a magnetic conpass ftom a geographical re- 
ference; measuring the depth of the spaced positions from a known reference 
point; assuming a value and direction from the earth's magnetic field at the 
positiims; and determining by machine calculations the distance and direction 
to the cased borehole using said measurements and assumed value for the 
earth's field and location for the cased borehole, said calculations using a 
least squares fit of the measured values to calculated values in the formula 
for the magnetic induction surrounding a long cylinder, 

2. The method of claim 1 wherein the radial co>q>onent of the field 

contributed by the casing is given by the formula: 

B - > a >^i{2niro/L) K, (2nirp/L) 

^ n - 1 " K.(Wl) f2nirZ/L) 4 /^^^nvti/l) PnirZ/L) 

and the axial component of the field contributed by the casing is given by 
the formula: 

B, . 'ZZ a„ ^0^^"^^^^ sin (2nir2/L) - b ^t2nirp/L ) cos (2niTZ/L) 
n»l " K^(2nira/LJ " K^(2niia/L) 

where a is the outer radius of the casing, Kq and modified Bessel func- 
tions of the second kind* a^ and b^ constants, L the wave length of the 
periodic magnetisation, p the distance from the axis of the cased hole to 
the point of measurement and Z the distance parallel to the axis. 



5. Hie method of claim 2 wherein L is the average length of one joint 

of casing* 
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4. The nethod of claim 2 whoroin L is larger than the interval over 
which data is being fit. 

5. The method of claim 2 wherein the measured data is substituted in 
the formula and the least squares fit of the data calculated to the measured 
data obtained, the calculated results then being used to determine the 
calculated deflection of a magnetic coopass, 

6. The method of claim 5 wherein the actual deflection of the magnetic 
coB5>ass is con^jared with the calcuUted deflection to determine the location 
of the cased borehole. 
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